Relic gravitons as the observable for Loop Quantum Cosmology 

Jakub Mielczarelfl 

Astronomical Observatory, Jagiellonian University, 30-244 Krakow, ul. Orla 111, Poland 

Marek SzydlowskQ 

Astronomical Observatory, Jagiellonian University, 30-244 Krakow, ul. Orla 111, Poland and 
Marc Kac Complex Systems Research Centre, Jagiellonian University, ul. Reymonta 4, 30-059 Krakow, Poland 

In this paper we investigate tensor modes of perturbations in the universe governed by Loop 
Quantum Cosmology. We derive the equation for tensor modes and investigate numerically effects 
of quantum corrections. This investigation reveals that the region of super-adiabatic amplification of 
tensor modes is smaller in comparison with the classical case. Neglecting quantum corrections to the 
equation for tensor modes and holding underlying loop dynamics we study analytically the creation 
of gravitons. We calculate the power spectrum of tensor perturbations during the super-inflationary 
phase induced by Loop Quantum Gravity. The main result obtained is the spectrum of gravitons, 
produced in the transition from the quantum to classical regime of the Universe. Obtained spectrum 
is characterized by a hard branch. The numerical investigation shows the strong dependence of ^ max 
frequency with respect to quantum numbers. The parameter f2 gw (t'max) approaches even to ~ 1CP 1 
for highest frequencies. We compare our results with recent constraints from the Laser Interferometer 
Gravitational-wave Observatory (LIGO) and find that it is possible to test the quantum effects in 
the early Universe. 



I. INTRODUCTION 



Loop Quantum Gravity (LQG) introduces strong modifications to the standard description of the early universe. 
The main difference with the classical approach is the avoidance of an initial singularity This effect leads to 
the bouncing solution on the semi-classical level [3, 3 ■ Another interesting property is the occurrence of the super- 
inflationary phase induced by quantum effects This phase is in fact not long enough to explain the observed 
flatness of the Universe, but after this phase the Universe has proper initial conditions to start the standard slow-roll 
inflation. In this scenario a inflaton field firstly climbs up the potential hill and then stops before a slow-roll phase, 
producing the running of the spectral index [5[. The production of scalar perturbations during super- inflationary 
phase is investigated in the papers [f| 0, @] • 

In this paper we consider the transition from the quantum to classical universe through the super-inflationary 
phase. Because during this transition the main contribution to the energy of the inflaton field comes from the kinetic 
part, in calculations, we neglect the contribution from the potential energy. It is worthwhile to note here that the 
super-inflationary phase inducted by LQG is a generic property and does not depend on a kind of the field which fills 
the universe. For analytical considerations we solve dynamical equations in the semi-classical and classical regimes 
and then we match them. It is done for the value of the scale factor 

a = a* = y^^ pl W 

where j is a half-integer quantization parameter, lp\ is a Planck length and 7 is the Barbero-Immirzi parameter. The 
latter parameter 

m 2 

comes from calculations of black-holes entropy Q. Below the value a* non-perturbative modifications become im- 
portant. We derive the equation for tensor modes in the LQG scenario. We investigate numerically effects of loop 
corrections. Neglecting quantum corrections to equation for tensor modes and holding underlying loop dynamics 
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we study analytically creation of gravitons. We calculate the spectrum of tensor perturbations during the super- 
inflationary phase and the density of gravitons produced during the transition from the semi-classical to classical 
universe. Such gravitons give contributions to the stochastic background of gravitational waves. Nowadays the detec- 
tors like LIGO aim at the detection of these stochastic gravitational waves [ll[ . Usually to describe the spectrum 
of gravitational waves the parameter 



f2 gw (^) 



p c dv 



(3) 



is introduced. Here p c is the current critical density, /9gw is the density of gravitational waves and v is the physical 
frequency measured today. The recent LIGO constraint for this parameter is f2 gw < 6.5 • 10 -5 [12j. We also calculate 
the value of the function fl gw (v) in the model and compare it with LIGO constraints. 



II. THE SEMI-CLASSICAL DYNAMICS 



Loop Quantum Gravity introduce strong modifications to the dynamical equations in the semi-classical regime. 
These modifications come from the expression for the density operator [131 ] 



d j (a)=D(q)- 



(4) 



where q is defined as follow 



a 

1 = — 

at 



and for the semi-classical universe (Zpi < a <C a*) the quantum correction factor has a form [D. 



D(q) = q^ 



1 



21 \l + 2 



[(q + 1)'+ 2 - \q - 1|'+ 2 ] - j^j [(q + l) i+1 - sgn(g - l)\q - 1|^] 



3/(2-2/) 



(5) 



(6) 



Here I is the ambiguous parameter of quantization constrained by < I < 1 [l5| . The Hamiltonian for the scalar field 
in the flat FRW universe has a form 
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H = ^dj(a)p^ + a 3 V{(t>) where = d } 1 (a)4> 



(7) 



This lead to the equation of motion of the form 



(8) 



The Friedmann and Raychaudhuri equations for the universe filled with a scalar field are respectively 



H z = 

d 
a 



8ttG 
8ttG 



2D 



+ v(4>) 



D 



D 



4HD 



From equations ([8]) and ([9]) we obtain the relation 



(9) 
(10) 



H = -4ttG— 1 — 

D \ 6HD 



(11) 



Due to quantum correction D in the region (lp\ < a -C a*), the expression in the bracket can be negative, leading to 
H > (super- inflation). If a ^> a* then fl«l leading to H < (deceleration). For a <C a* the approximation of 
expression ([5]) have a form 



D(q) 



3 v 3/(2-2/) , a N 3(2-0/(1-0 

TTT 



(12) 
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We use this approximation to calculate the dynamics in the semi-classical region. Now 

D _ 3(2-0 



HD l-l 



> 6 (13) 



leading to the phase of acceleration, see equation (fTTj) . Putting (fT3"|) into equation (fTTj) and combining with © we 
obtain the equation for the scale factor 



{a'f 



n 3 1 
2 + - 



2 1-/ 



= (14) 



where prime means the derivative in respect to the conformal time dr = dt/a. We assume here V((j>) = as it was 
mentioned in section I. The solution of (TTJJ is of the form 

a oc (-tT 2 5Tt. (15) 

To calculate the solution in the classical regime we take D = 1. In this case the equation for the scale factor have a 
form 

aa" + (a') 2 = 0. (16) 

Now we match two solutions, from two regions, at some To as follow 

ai(-To) = 02 (-to) , (17) 
a'i(-To) = 4(-to). (18) 

Where region 2 is classical and region 1 is semi-classical. The value of the chosen conformal time tq corresponds to 
the scale factor a* . After matching we obtain the solution of the form 

for t < -t , (19) 

1 for t > -t . (20) 

This solution is shown in Fig.Q]together with the numerical solution. The upper curve corresponds to the evolution of 
the scale factor o(t), while the bottom curve does to the first derivative of the scale factor in respect to the conformal 
time. The obtained solution is of course only an approximation of the real evolution, however it is sufficiently exact 
to be used in analytical calculations. We find the agreement with the numerically calculated evolution of the scale 
factor and the Hubble rate obtained by Tsujikawa et al. @. In the future investigations we use both numerical and 
approximate analytical solutions. 




III. EVOLUTION OF TENSOR MODES IN LOOP QUANTUM COSMOLOGY 

Tensor perturbations hij to the FRW metric we can express as 

ds 2 = o 2 (t) [-dr 2 + (Sij + hij)dx % dx j ] 



(21) 



with \hij \ <C 1. Using constraints h\ — ^ihj — we can see that tensor hij have only two independent components 
h\ = —h\ = h + and h\ = h\ = h x . These components correspond to two polarizations of gravitational waves. Since 
tensor modes of perturbation are not coupled to the scalar field source, we can obtain equations for them from the 
variation of the action 



S, 



(2) _ 



1 



64ttG 



d i xa 3 



1 



1 



32ttG 



d^xa 3 



1 / 



a' V 



1 



(22) 



For the detailed discussion of this issue see Ref. [16|, llTj . Both polarizations of gravitational waves are not coupled 
and can be treated separately. Introducing 



(23) 
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FIG. 1: The evolution of the scale factor a (upper curve) and a' (lower curve) in the conformal time (with to = 5, j = 100 
and I — 0.1). The dashed line presents the approximate solution and the solid line corresponds to the numerical solution. The 
initial conditions are fixed for a,. 



we can rewrite the action for tensor modes in the form 

4 2) = \ I rfv 



(24) 



Quantum corrections can be introduced now in the same way as in the scalar field case , [l4[ , leading to the equation 
of motion 



\ D a A 



(25) 



There are no other corrections from Loop Quantum Gravity which contribute to this equation. 

Introducing the new variable fj, — ah and changing the time for conformal time we can rewrite equation (|25[) to the 
form 



-DV 2 



a a D 



(26) 



Since the fluctuations considered have the quantum origin we must change the classical y, for the corresponding 
operator //. The field fx and conjugate momenta fr can be decomposed for the Fourier modes according to 



fi(x, t) 



d A k 



2 ( 2*)3/*y-"{^ + ^* 3 }' 



(27) 
(28) 



where the relation of commutation [jl(x, r), ^(x, r)] = iS^(x — y) is fulfilled. The equation for the Fourier modes is 
now 



D' 



(29) 



where 



„ r2 1 fa" a' D' 
M = — — 



D \ a a D 



(30) 
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is called the pump field. In the classical limit (D — 1) equation ([29]) assumes the known form 

H = °- ( 31 ) 

Because it is impossible to solve equation (|2"5|) analytically we must investigate the effect of quantum corrections 
numerically. Because D is always positive we can have amplifications of the tensor modes when k 2 < M 2 . The pump 
field function was shown in Fig. [2] In the first panel (left up) we draw pump field M 2 with neglected quantum correc- 
tions calculated numerically and with the use of solution (|20j) . What we see is that the numerically calculated pump 
field extends the region of super-adiabatic amplifications. In the next panel (left right) we compare the numerically 
calculated pump field with and without quantum corrections. We see that quantum corrections lower the region of 
amplification. The obtained value is however still larger than this obtained using the approximated analytical solution 



a 
a 




FIG. 2: Top left: pump field M 2 with neglected quantum corrections to the equation for tensor modes calculated numerically 
(circles) and with use of solutions (|20|l (line). Top right : pump field M 2 with quantum corrections to the equation for tensor 
modes (boxes) and without corrections (circles). Bottom left : Global behavior of the pump field M 2 with quantum corrections 
to equation for tensor modes. Bottom right : Evolution of the friction term D' / D in the equation for tensor modes. In all 
panels it is assumed tq — 5, j — 100 and I = 0.1. 



The friction term D' /D in equation (|29|) can be positive or negative leading to the amplification or to the damping. 
We see that the friction term D'/D is for the most of time positive leading to amplification and drastically peak to 
negative values in the neighborhood of a*. In fact this dependence strongly depends of a quantum number I. The 
increasing value of I the negative pick goes toward to positive values and becomes less sharp. 

We can now use the approximation D = 1 in equation ([29]) to calculate the spectrum of tensor perturbations during 
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the super-inflationary phase. The spectrum of tensor perturbations can be now expressed using correlation function 

,fi, H .,„. GAttG f d 3 k IA . .. its f dk^ ,,.sinfcr .„„. 
(0\h)(x,r)hl(y,r)\0) = —J J^f^)? ^ = J ~k ^~kr~ (32) 

where an Einstein convention of summation was used on the left side. For considered super-inflationary phase, using 
expression (|19[). we have solution 

K = ^V^HW(-kT) (33) 

where 

N = ^l e ^+i/2)/2 and (3 = 2^ (34) 

Normalization is found by correspondence to well normalized plane wave e~ lkT / 'y2k for high energetical modes 
|fcr| <C 1. Since for us interesting are super- horizontal modes we can use approximation 



HV{-kT)^- l -Y{v)(- K -L) . (35) 



Super-horizontal modes are these which firstly cross out the horizon and then evolve "frozen" in super-horizontal 
scales. Finally, in further epochs on universe, such a modes reenter horizon. These modes bring an information from 
earliest stages of the universe. Given modes cross the horizon when 

k ~ aH = £ (— Y , (36) 



t \a 

where we used definition of Hubble radius and equation of evolution ([19]) . Finally with use of definition ([32]) the 
spectrum at horizon crossing have a form 



where spectral index is equal 



and normalization constant is expressed as 



V T (k) = A 2 T k nT (37) 



UT = ~~dlnk~ = 2T7 (38) 



In this case tensor spectral index (|38[) is positive and tit S (0,3). To compare, from the standard slow-roll inflation 
tensor spectral index is tit = — 2e + 0(e 2 ). Unfortunately available data from CMB or from large scale structures 
observations are not sufficiently precise to determinate value of the tensor spectral index. 



IV. RELIC GRAVITONS FROM THE QUANTUM TO CLASSICAL UNIVERSE TRANSITION 

In the previous section we derived and investigated equation for tensor modes P^]). Using approximations we also 
calculated the spectrum of gravitons which cross the horizon during the super-inflationary phase. Now we want to 
calculate a number of gravitons which are produced during the transition from the quantum to classical regime of 
evolution. Before we start it, let us calculate the width of the band of produced gravitons. It can be directly taken from 
condition k 2 < M 2 . We use the classical approximation for the evolution of tensor modes for it preserves the physical 
picture of the graviton creation process as the previous numerical investigations indicate. In this approximation 
(D = 1) we have a maximum of frequency for tq, so with the use of equation (|19|) we have 



u = vW + i)- 



(40) 
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In fact, as it can be seen in Fig. [2 this frequency is generally higher. The corresponding maximal frequency for the 
present epoch is 



_. max = I ) = : ___ ____ i i . (41) 

27ra» \a to day/ 2ira* t Va to day/ 

To estimate this value we can approximate 

"today _ JP1_ = 1.4-10 32 K ^ iq32 

a* Tcmb 3.7 K 

where Tpi is the Planck temperature. The other way to estimate value atoday/a* is to use the Friedmann equation 
with radiation. This gives equation 



VW + 1) 1 ( a* 



Otoday \ H today HtodaytpiTQ Hj 

^r = v _ ^r = V — p — vi (43) 

where we use solution (|20[) to calculate iJ*. To obtain a numerical value we must know To. We use here the constraint 
for energy in the form l^l/mpj < 1 [l_| (the kinetic energy dominates over the contribution from the potential part 
as we mentioned before) for ai — ^/7_pi. Below the value of chosen at the space becomes discrete and the smooth 
dynamical equations cannot be used. The boundary for the kinetic energy is introduced to prevent energies beyond 
the Planck scale being produced. With use of this boundary conditions and the Friedmann equation ([9]) with solution 
(|2H|) we obtain the constraint for the conformal time To 

^—A-X—i) (j) - (44) 

As an exemple for the model with I = 0.1 and j = 100 we obtain To > 0.0014 and for the model with I = 3/4 and 
j = 100 we obtain To > 1.6 • 10~ 8 . We see that this boundary depends very strongly on the quantum numbers. 

______ _____ 1/2 

Combining equation (|4"Tj) with (|4"3"|) we see that f max oc t , so the boundary (|4"4")1 gives us also the upper constraint 
for a maximal value of frequency z^ max . For the model with I = 0.1 and j = 100 we have v max < 6.6 • 10 Hz and for 
the model with I = 3/4 and j = 100 we obtain f max < 2.8 • 10 24 Hz. Generally values of f max can be smaller than 
boundary values. For the further studies we choose the model with to = 0.1. So in this case the width of the band of 
relic gravitons considered is nowadays [0, 10 3 GHz] for I = 0.1. 

Fourier modes of (J27J) and ([__%)) for the super-inflationary evolution (TT9"|) can be written with the use of annihilation 
and creation operators as follow 

Afc( r ) = Hh{k,T) +al j .f*(k,r) for t<-t , (45) 
%( T ) = H9i{k,r) +a f t g*(k,T) for r < -r . (46) 



In this case the values of coefficients are 

\ 

hk 

gi(k,r) = - 

where 



A(fc,T) = ^V=kFHP(-kT) , (47) 



-HlU-^ + ^^H^i-kT) 



(48) 



A/". = y| e Wi/2)/2 and v = /3+l (49) 

Similarly, modes of (j!?T)) and for the classical evolution (T__TJ|) we can be written down as 

AfeW = 6fc/2(fc,T)+6 t _-/*(fc,T) for t>-t , (50) 

%(t) = b % g 2 {k,r) +P_^(k,r) for t > -t . (51) 
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Where the coefficients of decomposition are 



/a(*,r) 



92{k,r) 



with 




7V2 

TO 



2 + 1 



(fer + fcC) 2(1 - I) 



70 



TO 



1 



A 1 — 1 ( Tq+T 

2+1 \ to 



(fcr + fcC) exp (ikC,) , 



2 + 1 



- kr Q \ 



2 + 1 



TO 



t )h[ 2) (kT + kC) 



(52) 



exp {ikC) (53) 



2 + 1 



1/4 and ( = t - 



32-1 



4 1-/ 



(54) 



where is the Haenkel function of the second kind. 

The main idea of particles creation during transition comes from the Bogoliubov transformation 

b % = B + (k)a x + B_(kra^, 
bt = B + (kyal + B_(k)a_£ 



(55) 
(56) 



where from relations of commutation [fig, at] = 5^(k — p) and [6g>&t] = S^(k — p) we have \B + \ 2 — \B-\ 2 = 1. In 
the quantum phase we have d^\0i n ) = where \0i n ) is the vacuum state of this phase. In the final classical epoch, 
similarly bc\Q ou t) = what defines the new vacuum state \0 O ut)- But since we are in the Heisenberg description the 
true vacuum state in the classical phase is |0j n ) and thanks to the mixing from the Bogoliubov transformation (|55[) 
we have &g|0.j n ) = B^(k)*a ! g|0.j n ). So when B_(k) is the nonzero coefficient we have the production of particles 

(gravitons) in the final state. What we need now is to calculate coefficients of the Bogoliubov transformation B-{k) 
and B + (k) which can be written as 



B-(k) 
B+(k) 



/l(-T0)fl2(— n>) - gi(-To)/ 2 (-To) 

fi (-T )S2(-T ) - 52(- T o)/2(-T ) 

Ii(-tq)92{-tq) - gi(-r )/ 2 *(-T ) 

/2(-To)52(- T o) - ff2(-T )/ 2 *(-To) ' 



(57) 
(58) 



Since the total momentum of produced gravitons is conserved we can write the expression for the number of produced 
particles 



(59) 



As we can see, to calculate a number of gravitons we only need to know the coefficient B-(k). Now we can calculate 
the function O gw (z^) defined in equation ([3]). The essential energy density is from the relation 



dp gw — 2 ■ Hlu 



A-KLU 2 d(jJ 

{2-Kcf 



(60) 



where factor 2 comes from two polarizations of gravitational waves. With the use of relation (|4"Tj) we finally obtain 
the equation 



O gw ( I y) = 3.7 • l(T 49 /io V n (V/3(/3 + l) — 



(61) 



where ho is the normalized Hubble rate ho = -ffo/100 km -1 s Mpc. We compute this function and show it in the 
logarithmic plot with I = 0.1 and 1 = 3/4 (Fig. [3]). This spectrum is characterized by a hard branch with the maximum 
for ~ 10 12 Hz for / = 0.1 and ~ 10 11 Hz for I = 3/4. In this limit f2 gw approaches respectively to ~ 10 _1 and ~ 10 -5 . 

When the high energy region is shown only, the dependence S! gw (i/) on the quantum parameter I = 0.01, 0.1, 3/4 is 
exhibited (Fig. g]). 
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FIG. 4: Function fi gw (i/) for / = 0.01, 0.1, 3/4 (from top to bottom), to = 0.1 and j = 100. Frequency scale in Hertz. 



As we mentioned in section I, recent constraints from LIGO are O gw < 6.5 • 10 -5 [lj]. The LIGO observations are 
however concentrated in the region of ~ 10 2 Hz. From Loop Quantum Cosmology we have in this region f2 gw ~ 10 -14 
(for I = 0.1), what is extremely below the observational sensitivity. The numerical values obtained by us contain 
estimations of the time of transition to classical universe. The used value should be somehow proper to the order of 
magnitude. So we expect also similar deviations of fl gw (is). 

The spectrum obtained here is not a distinct feature of Loop Quantum Cosmology. As it was shown by Giovannini 
(T^ | a similar high energy branch was obtained in the quintessential inflationary model. Calculations based on String 
Cosmology lead also to similar results [2(|. To compare, for the standard inflationary models the spectrum is flat. 



V. SUMMARY 



Loop Quantum Cosmology has received much attention in the theoretical astrophysics. But what was lacked so far 
was empirical consideration of this theory. Bojowald indicated the quantum effects are negligible small at the present 
epoch but they can potentially tested [2l|. Along Bojowald's lines we showed that gravitational waves can be the real 
observable for testing Loop Quantum Gravity effects. 

We have considered the transition from the semi-classical to classical universe described by Loop Quantum Cos- 
mology. In the analytical approximation we obtained the tensor energy spectrum of the relic gravitons from the 
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super-inflationary phase. The analytical model takes into consideration the corrections to dynamical evolution only. 
While taking corrections to the equation for the tensor modes this equation cannot be solved analytically, so it is 
only possible to consider it numerically. The numerical investigation of the equation for tensor modes gave us that 
lower !/ ma j is admissible when the loop quantum effects are incorporated. The loop quantum gravity effects product 
additional damping during the production of gravitons. This is a challenge for future investigation - the full numerical 
analysis of this model. 

When we considered the production of gravitons during the transition phase the spectrum of these gravitons 
is characterized by the hard branch. The corresponding value of the parameter f2 gw , in its maximum, is f2 gw ~ 
10~ 7 . . . 10 _1 , depending on the value of the parameter of quantization I. In the region of the LIGO highest sensitivity 
we obtained the very small value of the parameter f2 gw , namely ~ 10~ 14 for I = 0.1 and ~ 10~ 28 for I = 3/4. As we 
mentioned, the similar hard branch is also a feature of quintessential inflationary and String Cosmology models. This 
work gives the motivation to search for high energetic gravitational waves. 
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